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1.0 SUMMARY 

A new exoatmospheric, powered explicit guidance (PEG) thrust integral 
formulation and a simple method of implementation is presented in 
this note. The new thrust integral formulation is significantly 
simpler than that currently used in PEG. Preliminary estimates 
indicate a computer storage savings of 220 words, which is approxi- 
mately 10 percent of the current PEG ascent program. Alternate 
methods of implementation that could produce even more savings are 
noted . 

The method of implementing the thrust integrals derived herein does 
not represent a departure from the current PEG implementation approach. 
Consequently, required simulation verification of the new equations is 
minimum. Essentially, one equation replaces the higher order thrust 
integrals of the current PEG equations, and this one equation can be 
verified analytically. The current PEG higher order integrals (higher 
than first order) formulation results from approximating the sine and 
cosine functions as series expansions. These approximations are 
eliminated in the formulation presented in this note. Therefore, the 
thrust integral formulation proposed herein is more accurate than the 
current formulation, when large steering angles are involved. A 
simple function or factor (f^) is developed such that the higher order 
integrals are a product of the first order integrals and the quantity 
( 1 -f -j ) . The equation for the factor f-| is. verified analytically. The 
first order integrals of the current PEG equations are simplified and 
maintained, and therefore, require no ''erification. 



- 1 - 



f 

\ 


1 


The proposed equations and method of implementation offer the flexi- 
bility of applying a different guidance law by changing the value of 
one scalar parameter. This additional flexibility, producing more 
accuracy and optimality, is gained at significantly less cost, and 
is implemented in such a manner as to have no impact on other GN&C 
programs, e.g., the G&C steering interface routine 

The new thrust integral formulation was tested and verified in simu- 
lations of Baseline Reference Missions 1 and 3A. Each simulation of 
the new integrals yielded virtually the same trajectr-^y and performance 

i 

as the current integrals, assuming the same guidance law (e.g., linear 1 

angle steering). However, the method of implementation proposed in this 

note allows high performance missions such as Mission 3A, to fly linear 

tangent steering. The new thrust integral formulation coupled with the 

proposed implementation method yielded slightly better performance 

than the current PEG equations for Mission 3A (i.e., +70 pounds MECO 

weight). j 

The new thrust integral formulation presented in this note reduces 
flight computer memory and processing loads for PEG guidance and is 
readily implemented within existing PEG. It is recommended that this 
thrust integral formulation be incorporated into the Shuttle powered 
flight guidance software requirements. 


- 2 - 



2.0 INTRODUCTION 

The Mission Planning and Analysis Division, supported by MDTSCO, 
is engaged in the development and verification of guidance soft- 
ware requirements for implementation in the Shuttle onboard GN&C 
computers. This activity is being conducted as established in 
Track Task Agreements between Rockwell International and the 
Shuttle Program Office. 

Development of the new powered explicit guidance thrust integrals, 
described herein, was initiated by request of the Shuttle Powered 
Flight Guidance Working Group, chaired by Aldo Bordano, FM7, in 
support of the guidance scrub activities started in April 1976. 
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3.0 DISCUSSION AND ANALYSIS 

This section is divided into six subsections preceded by a discus- 
sion of the exoatmopsheric powered explicit guidance problem and 
the relationship of thrust integrals to the solution of this 
problem. The new thrust integral formulation is developed in the 
six subsections. Subsection 3.1 provides a slightly simplified 
development of the basic first order thrust integrals currently 
used to solve the powered explicit guidance problem. Since the 
total thrust integrals are integrals of the thrust acceleration 
vector, the equation defining the optimum unit thrust vector is 
derived in Subsection 3.2. The approach used in generating a new 
thrust integral formulation is explained in Subsection 3.3. Partial 
integrals, i.e., assuming that the time about which the guidance 
solution is expanded is one half of the time-t'o-go (K = Tgg/2), are 
developed in Subsection 3.4. Total thrust integrals, where K = T^q/2 + 
AK (AK is a small perturbation term), are developed in Subsection 3.5. 
The proposed method of implementing the new thrust integrals, including 
flow charts, is presented in Subsection 3.6. 

The basic rocket vehicle vacuum flight equation of motion is of the 
form 

V = ^ If + where (1 ) 
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* 

V * Total acceleration vector, 

F = Engine Thrust, 

m = Vehicle mass, 

F/m * Thrust acceleration, 

= Unit vector in desired thrust direction, and 
G = Gravitational acceleration vector. 

The typical powered flight guidance problem is to determine real 
time values of i^ (the unit thrust vector) and the rate of change 
of i^ that steer the vehicle to desired targe;, conditions while 
expending minimum fuel, i.e., the fuel optimum path is determined 
between current state and target state. 

Explicit guidance provides a closed-form predictor/corrector solu- 
tion to the two point boundary-value problem discussed above. 
Closed-form integrals of the thrust acceleration vector, (F/m)i^, 
are involved. The thrust integrals are not exact, in general. 
However, thrust integrals are developed in the subsequent sections 
that produce more accuracy than is necessary in practice, since the 
integrals converge to exactness as the remaining burn time, Tqq, 
approaches the value of zero. The thrust integrals developed in this 
note are exact for a constant acceleration burn and virtually exact 
for a low thrust burn. The integrals are simple and provide a 
qualitative understanding of the exoatmospheric trajectory optimiza- 
tion and powered flight guidance problem. No rigorous mathematical 
proof is attempted in this note, i.e., a simple and practical 


- 5 - 







l 1 1 1% mrniku 


t If i>*mfciiniirtii 


engineering approach is taken for constructing a practical solution. 

In Section 3.1 it is shown that by making small angle assumptions, 
the unit thrust vector can be approximated in the form * a. + bt, 
where t is time and a^ and b are constant vectors, and the thrust 
acceleration vector (Vy) assumes the form Vy = — (a^ + bt). 

Integration of the above equation involves integrals of the quantities 

£ and £t. These are conventional first order thrust integrals, 
mm 

Multi-stage equations for these basic thrust integrals are developed 
in Section 3.1 . 

In Section 3.2 it is shown that the fuel optimum acceleration vector 
is approximated in the form 

lj = m Ci cos + (b/u>) sin uit], 

where a^ and are constant vectors and to is a constant angular rate. 

In practice, the independent variable time, t, in the above equations 
is redefined as z = t - K, where K is the time about which the 
solution is expanded and the value of t ranges from zero to time- 
to-go, T Gt) , i.e., the integration limits are from zero to Tqq. The 

* 

value of K is a function of the thrust integrals. Assuming constant 
acceleration, K = Tqq/ 2 and in any case K*=Tqq/ 2. The explicit 
guidance problem reduces to determining values of K, a, and b that 
satisfy the desired target renditions (in the following sections, a 
and b_ are redefined as \ and X). When dealing with the sine and 
cosine functions, it is convenient to define K as K = Tqq/ 2 + AK and 
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solve for a value of AK, 1.e.» the following expansions are employed: 


cos (x - 6 ) ■ cosx cos6 + sinx sin6 and 
sin (x - 6) = sinx cos6 - cosx sin6, where 
xs w(t - Tqq/ 2) and 6 = wAK 

* 

Integration of the above equation for Vq. involves integrals of the 

quantities £ cos ut and £ S1 ‘^ . u>t . Assuming constant thrust and mass 

flow rate, closed form integrals of these quantities do not exist. 

2 2 

In the current PEG formulation, it is assumed that cos wt = 1 -ut/2 

and that (sin wt)/w = t. Therefore, the first order integrals, as 

F 2 

developed in Section 3.1, as well as second order integrals of -t 
are utilized in the current PEG equations. The approach taken in this 
note eliminates these second order integrals. 

F F 

Closed form integrals of - i^ exist if - is expanded as a polynomial. 

The integrals involve integration of quantities of the form x n cos x and 
x sin x, which have closed-form solutions. It is seen that the solution 
is exact assuming constant acceleration or linear acceleration, as is 
nearly the case for low thrust OMS burns. And in any case, the solution 
converges to exactness as Tqq approaches zero. 

Using mean value considerations, it will be shown in Section 3.5 that 
the integrals of ~ cos wt can be adequately represented by assuming a 
constant, mean acceleration. However, this is not the case when dealing 
with the quantity — It is shown that the first integral of this 

quantity is adequately expressed by assuming linear accelt 'ation. It 


*Where n is the order of the polynomial. 
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Is illustrated that the second integral of £. sJ_n_cot_ re q U ^ res an 
order greater than first*. However, this integral can be 
adequately expressed in terms of the first order integrals developed 
in Section 3.1. Consequently the integrals developed in the follow- 
ing sections assume a constant, mean acceleration or a linear 
acceleration profile. A second order acceleration profile is dis- 
cussed for illustrative purposes. 


*i.e., if j-'j' is expanded as a polynomial. 
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3.1 Basic Thrust Integrals 

The basic first order integrals used for solving the powered explicit 
guidance problem are developed below. The integrals are the same as 
in the current PEG equations except that they are somewhat simplified. 

For the shuttle vehicle, either constant thrust (and mass flow rate) 
or constant acceleration is assumed. If constant thrust is assumed, 
the acceleration magnitude as a function of time is 

% 

a (t) = F/(m - mt), where 
0 

t = time, 

F = constant thrust, 

m Q = initial mass, 

m = constant mass flow rate. 

The above equation can be written in the form 

a(t) = (F/m)/[(m o /m) - t] or 

a(t) = V ex /(t - t), where 

F/m = V ex (constant exhaust velocity), and 

m 0 /m = x. 

It is convenient to define the parameter t as 

• • 

t = m o /m = (F/m)/(F/m o ) or 

T ■ W for 

-9- 


•a 



X 


Vg X one of the most constant parameters in the propulsion system 
(it is considered constant in the guidance problem) and the current 
acceleration, a Q , is measured by the navigation system, whereas m Q 
and m are not so easily determined. 

The basic, constant thrust integrals are presented here for the i th 
stage of an n stage vehicle. 


T Bi 


'Bi t 

L i ■ / a i dt - and s i 'LL a.dsdt, where 


o o 


a i = Wh - l > and 


= Burn time of 1 th stage. 

Employing integration by parts, the integrals are easily determined 
as 

L i = V exi ln M T i " V^exi ln ^ " T 81 /t 1> and 


^ = -L^t. - T b1 ) + V ex .T 


Bi 


( 2 ) 

(3) 


During the SSME constant acceleration burn 


a.. = a^ (constant acceleration limit), 


and the thrust integrals are simply 


L. = a L T B . and 


Si = .5L.T B .. 


(2-a) 

(3-a) 
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In equation 1, it can be assumed that if ** X. + X_ (t - K)*, where 


K *=>1^/2 (to be determined). 

viJ 

Tqo = Time-to-go, or remaining burn time. 

X^ = constant unit vector when t = K. 

* 

X^ = constant rate vector normal to X., i.e., 
X • X = 0, 

X • X = 1. 


It is assumed that XK is small compared to unity (small angle 

approximation). Further assume that burn times (T^.) of all stages 

are known, and that K, _X and \ have known values. Now, the thrust 
% 

component of V, in equation 1, can be integrated to produce velocity 
and position changes due to thrust (V^ and Rj), i.e., 

V T = “ lk + i (t - K)], 

V, = f G0 V T dt, and 
J o 

-t = f G ° f* i T dsdt - 

J o 'o 

The integrals of ^ (t - K) are presented in this section, in addition 
to the integrals of j~. These are conventional first order integrals 
used in the current PEG equations, and will bt employed in the 
simplified thrust integrals presented in this note. 


*This is known as linear tangent guidance (LTG). 
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Performing the above integration yields equations of the form 


r f 


Vj - LX + X ( J - LK) , and 
Rj = SX + X. (Q - SK), where 


(6) 


(7) 


L and S are as defined by equations 2, 3, 4, and 5 and J and Q are 
developed in the following paragraph. 

Using the same nomenclature as for and S^, 


T Bi 

J.j = / D a . ( t + t Q .j) dt and 
o 


Tp j t 

Q i = / / a i (s + t .) dsdt, 

o o 


where t . is the initial time of each stage (called t „ in the 
oi 3 ' goa.. 


current PEG equation) i.e., 
f Bj 


i-1 

tni = E T Ri ’ 


j=k 


e.g. 


V °* T o2 = T B1* T o3 = T B1 + T B2* etc ' 


Employing integration by parts, the constant thrust integrals are 
easily determined as 


J i ■ L i T i ' V exi T B1 * L i ‘oi- and 


■ s i '1 - 1/* V oxi T2 Bi + S i ‘or 
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From inspection of equation 3, it follows that 


J i 3 L 1 T Bi - S i 4 L 1 ‘oi 


or 


J 1 ' L 1 ‘goi ’ V where ‘goi s ‘oi + T Bi ' 
The equation for is written as 

Qi ■ S. (t. + t Qi ) 1/2 V exi -T| i . 


For a constant acceleration phase, a^ * a^ (constant accelera- 
tion limit), and the thrust integrals are 


( 8 ) 

(9) 


J i = S i + L i‘oi a " d 


Q,=l/3 S lTe1 + S ( t oj H S, (TB t /3 + t o1 ). 


where 


L i = a L T Bi and S i “ • 5 L i T Bi 


Equation 10 can be written as 


d 1 ■ S T B1 - - 5 L i T Bi + L i‘oi or 


J. = L.t . - S. 
l i goi i 


where T ^ is defined by equation 8. 


(10) 

(ID 


( 12 ) 
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Inspection of equations 8 and 12 shows that the expression for 
is the same for constant acceleration as for constant thrust. 
Summing over all stages (i = k, n), the multi-stage integrals are 


n 


E J i 

and 

(13) 

i=k 



n 

i-1 

(14) 

E Q, 

i=k 

+ where J of E Jj • 


It is easily verified that 

0 - E <4Vl - V = LT go - S, where 
1=k 

n 

^GO = £ ^Bi (Time-to-go) , and L and S are multi-stage 
i»k 

integrals as defined in equations 4 and 5, i.e., 

J = LT gq - S. (15) 

In practice, K, X^ and X_ are not known until after the thrust 
integrals are formed. It will be seen later that K is chosen such 
that the coefficient of X in equation 6 vanishes, i.e., 

0 - LK = 0 or 

K = J/L. 

Using this equation and equation 15 

K=T go -S/L. (16) 
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The first order integrals are summarized as follows: 
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i 


CONSTANT THRUST 


L i = " V exi ln ‘ WV 


S i = " L i ^ T i " T Bi* + V exi T Bi 


Q. = S. (t. + t 0 .) - J v exi T ! 


2 

Bi 


CONSTANT ACCELERATION 


L, = 


a L T Bi 


S i 


*i 


2 L i T Bi 


S i (f T Bi + t oi> 


CONSTANT THRUST OR ACCELERATION 


J. = L.t . - S. 
l i goi i 


( 2 ) 

(3) 

(9) 


(2-a) 

(3-a) 

(ID 


( 12 ) 


The respective multi-stage equations for L, S, J, and Q are 
equations 4, 5, 13 and 14. The above equations are somewhat 
more simple than the current PEG thrust integral equations in 
that one equation for applies to constant thrust or constant 
acceleration. 


i< i 
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3.2 Fuel Optimum Unit Thrust Vector* 

In this subsection, a well known method of the calculus of variations 
(the Euler-Lagrange method of multipliers) is employed to derive an 
equation defining the optimum unit thrust vector time history. 

Assuming a spherical earth, it can be shown by employing the Euler- 
Lagrange method of multipliers that the fuel optimum thrust vector 
time history is defined by the differential equation 

X - -w 2 X + 3 (/(X/u^) u^ , where (17) 

Unit (X) (unit thrust vector), and 
w 2 = p/R 3 U> is circular orbit rate at radius R = |RJ ) 
p = gravitational constant 
£ = Radius vector 
u^ = unit (R) 

The following definitions are made: 

u_ e X_ (t) i.e., X as a function of time. 

z = t-K, where K is the constant reference time. 

u^ (K) = X, constant reference unit vector at t * K, or z = 0. 

u (K) = X, constant reference value at t = K. 

*This is sometimes called the prioier vector. 

-17- 
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The following simplifying assumptions and definitions are made: 


X • i - 0 (i.e., X is normal to X). 


w 2 is constant. 

Up a unit (fl + Rq) a constant mean value of Up. 

Rq = Desired, terminal radius vector. 

6 s sin" 1 (X. • Up), constant mean value of the angle between 
the thrust vector and the horizontal. (It will be seen later 
that X is approximately the mean value of u. = 2 lU).) 

With the above definition and assumptions, equation 17 can be 
written as 


V = -w 2 u[ + 3 id 2 sin9 cos wz Up, where 3w 2 sine Up 
has a constant, mean value. The cos wz term is added to 
compensate for the changing direction of that term when it 
is not assumed constant. Also, X^ cos wz is the dominant part 
of the homogeneous solution to equation 18. 

The solution to equation 18 is 

_u = cos wz + (X/w)sin wz+^(sine)wz sin wz Up 
and 

• • 

U[ = -wX_ sin wz + X_ cos wz 

+£• sin G (w sin wz + w 2 z cos wz) Up 

Differentiation of equation 20 verifies the solution, u^ is 

sometimes called the primer vector and u^ the primer rate 

vector. 


(18) 


(19) 

( 20 ) 
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The term involving sin 0 is usually small compared to the other 
terms, and is ignored in the following sections. For the shuttle, 
second stage boost to orbit, this term would increase the initial 
pitch angle by about 2 degrees, which would have very little 
effect on the burnout weight. 

Considering the integral 

AV = / GG F_ unit (u) dt, it is seen that 
o m 

AV (fuel expended) is minimum if sin 0=0. 

This implies that a maneuver involving no position constraints is 
optimum if the primer vector is normal to the radius vector 
at the midpoint of the maneuver, since is the unit radius vector 
at the midpoint of the maneuver, and it will be seen later that 
the reference time, K, has the value of approximately Tqq/ 2, there- 
fore, X. is approximately the value of the primer vector at the midpoint. 

Inspection of equation 20 shows that the optimum primer rate can not 
have a constant value unless sin 0=0. Assuming that the final 
optimum magnitude of u[, u^, is known, equation 2C can be used to 
produce a value of X = |Xj to satisfy this final value, i.e., 

Mf • = u 2 ^ (Substituting z = T^q - K in equation 20). 

This produces a quadratic equation from which X s determined. 

The equation can be made very simple by replacing u^ with 
• • 

unit (X) [u^ • unit (X.)], i.e., only considering the component of this 
term along the X vector, and the solution is of the form 

X « Vf (u77"u) - D. 
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This could be implemented in an approximate manner by assuming a 
constant, average primer rate of 


X = ylH Uf, w) - D/2, 

and recomputing this parameter each guidance pass. 

• • n • 

If u^ * tot then f (u^, to) = to and X = to (1 - D/to), where the 

* 

term D contains w. The term D has the sign of X. • u^. The term 

is time varying and vanishes when Tqq = 0. Also, if sin 0=0, 

D = 0. The above equation can be implemented in an approximate 
# 

manner as X = u (1 - D/2u), where w = WVpRp 2 . 

If the term involving sin 9 is neglected, equations 19 and 20 
become 

u^ = X cos toz + (X/ to) sin ooz and (21) 

u. = -uA sin toz + X^ cos wz, where (22) 

z = t-K, 
u 2 = p/R 3 , 

X • X^ = 1 , and ' 

X • X = 0. 


*i.e.» u =yfi/R 3 - p/V D R D 2 , where Vp and Rp are desired velocity and 
radius magnitudes. 
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Inspection of equation 21 yields interesting qualitative 
information: 

a) If u) = X, then u • u = 1 , i.e., u is a unit vector and the 
steering angle is linear with time. 

b) For a flat earth assumption (i.e., w approaches zero in equation 

21), equation 21 becomes 
u = X + X (t-K), 

which is the well known linear tangent steering law. If X is 

large compared to w, the above conclusion is approximately 
true. 

c) If there is a position constraint and X is not related to w, 
then the optimum steering is not necessarily linear angle or 
linear tangent. 

In the following analysis, for simplicity in understanding the 
concept, it will be assumed that w = X. Later, it will be seen 
that for Shuttle second stage burn this is not necessary unless 
throttling is involved (i.e., Return to Launch Site abort (RTLS) and 
Baseline Reference Mission 3B). For example, experience has shown 
that linear tangent steering is somewhat more optimum than linear 
angle steering in cases where X >w and throttling is not involved. 

This is especially true in Baseline Reference Mission 3A where the 
guidance assumes that an ngine is going out at the abort mode 
boundary. Linear tangent steering does not loft the trajectory as much 
as linear angle, and this is more optimum (if an engine does not 
go out). However, in Mission 1 where there is no pseudo engine out 
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the performance difference between linear steering and linear 
tangent steering is much less (*30 lbs. difference in weight 
at MECO). Saturn guidance used linear steering, which produced 
an orbital ^sertion weight loss of 100 pounds out of 250,000 
pounds compared to the theoretically optimum solution. 

• • • 

If it is assumed that w * X (where X = |x|), equation 21 be 
comes 

if = X cos X (t - K) + (yx) sin X (t - K), (23) 

and if * if = 1. i.e., if Is a unit vector. 

Equation 23 is the assumed form of the optimum unit thrust vector 
that will be used in developing the total thrust integrals. 

Since equation 23 produces a very nearly constant value of X 
throughout the flight, it can be assumed in developing the total 
thrust integrals (in the following analysis), that the converged 
value of X from the last guidance pass is used in defining the 
total integrals for the current pass. This is an important assump- 
tion since X for the current pass is not known, if a position 
constraint is involved, until the value of the vector, is 
determined, and the value of X. is a function of the total 
integrals. 
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3.3 Approach to Development of Total Thrust Integrals 

The approach used in generating a n*w thrust integral formulation 

is explained in this subsection. 

• 

Assuming that values of Tqq and X are available, and that the total 
thrust integrals (here defined as Ly, Jy, $y, and Qy) are known, 
the thrust t .locity and position assume the form, 


Vy = LyX + Jy X 

* a(t) i f dt, and 

0 

(24) 

Ry = SyX + QyX 

T rn t 

= / bU / a(s)i_ydsdt. 
0 0 

(25) 


The total equation of motion is 

V = Vy + G (where Vy = a(t)iy and G is the gravity vector). 

The following integrals are now considered: 

T G0 • T G0 . 

/ V dt =/ (V T + G) dt and 
o’ o 

/ G0 / Vdsdt = / G0 / (Vy + G) dsdt. 

0 0 0 0 

The final conditions are; when t = Tg Q , R = R^ (desired radius) 

and V = Vyj (desired velocity), and when t = 0, Bo E H ancl 1 q = I* 


Performing the above integration results in 


Bo - <£ ♦ vt go ) 

T 


• Br + VaV and 

■ B T ♦ Bgrav. whore 


V 

~grav 


/ G0 Gdt and ^ e/ go / Gdsdt. 
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For the development in this note, it is only necessary to assume that 
these gravity integrals exist. It is currently planned to make use 
of the numerical (average G) integration package (already in the 
navigation system) in order to extrapolate to the final state and 
obtain gravity effects over the trajectory, i.e., Vg rav and £g rav . 


Rearranging the above equations and using equations 24 and 25 for 
Vj and Ry results in 

LyX_ + JyX = and 

SyX + QyX = Rg N , where 

»fi N = - (v + lg r av > (velodty-to-so) and 

Bsn = Bo - (& + V T G0 + Vav> 

It can be shown that, from the standpoint of optimality, it is desir- 
able to expand about the vector and have no turning rate control 
in that direction (i.e., X. • = 0). Therefore, it is desirable 

to determine a value of K such that Jy = 0. 

Referring to equation 26, it is seen that if Jy = 0, then 
X - Unit (V™) and X • = 0, since X • X =0. 



i 

I 


Performing the vector dot product of X and equation 27 It follows 


that Sy * X • Rg N , and equation 27 is used to solve for X as 
X * / ^GN - x\, which Implies that R™ in equation 27 is 

7 

replaced with Syj ^GN Assuming that all necessary values are 

TO 

known (including K)> this completes the guidance solution for X^ 
and X. 


The reference time K can be defined as K = Tgg/2 + AK * The t,irust 

acceleration vector is written as 

• » • • * 

Vy = a[X cos Xz + (X. /X) sin Xz] or 

• » • * * 

Vy = (a cos Xz) X^ + [(a/X) sin Xz] A_, where z = t-T G Q/2 -AK, 
from which it follows (referring to equations 24 and 25) 
that T 


'60 


T G0 t 


Ly = / a cos Xz dt, Sy = / / a cos Xzdsdt 

o oo 


l G0 


T G0 1 


Jy = / (a/X) sin Xzdt, Qy = / / (a/X) sin Xzdsdt. 

^ o o 


t 


i 

Li I 

fe 1 


i 


If small angle approximations are made in the above equation 

for Vy, then 

• » 

Vy = a(X + xz) 

and the thrust integrals are the first order integrals of Section 
3.1, i.e., Ly = L, Sy - S, Oy * J-LK, and Qy « Q-SK. 
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The purpose of the next two subsections (3.4 and 3.5) is to determine 
K (or AK) such that Jy » 0, and develop functions (F^, an< * ^3) 
such that 



a cos Xzdt = F, / a dt * F.L, 

1 0 

t • T r t 

/ . a cos Xzdsdt = F « / adsdt * F^S, and 

0 J o 0 J 

/ (a/X) sin Xzdsdl = F ? Z 60 / azdsdt = F ? (Q-SK) , 
0 c o 0 c 


where L, S, and Q are the first order integrals of Section 3.1. 

F^, Fy,* an< * ^3 are simple function of X, Tqq, and AK. 

Therefore, the higher order thrust integrals of the current PEG 
equations are eliminated by introduction of the factors Fj, F^, and F^. 


In Section 3.4, partial thrust integral s* are developed assuming 
AK = 0, and in Section 3.5, total thrust integrals are readily 
derived from the partial integrals where 

x = X (t - Tg Q /2), 6 = XAK, and 

cos (x-6) = cos x cos 6 + sin x sin 6, and 

sin (x-<$) = sin x cos 6 - cos x sin 6. 

It is important to remember when reading Section 3.4 that, for conven- 
ience, K is defined as K = Tqq/2. In Section 3.5, K is defined 
as K = Tqq/ 2 + AK and AK is chosen such that 
T fn • 

J T -f (a/X)(sin x cos 6- cos x sinS)dt = 0. 

' 0 


♦Partial thrust integrals are here arbitrarially defined for convenience 
as integrals resulting from the assumption that AK = 0, or K - Tp n /2. 
These integrals are partial in that, in general, AK / 0, and u 
additional terms are involved. Total integrals include the additional 
terms when it is assumed that AK f 0. It will be seen that total 
integrals are simple functions of the partial integrals and AK. 



3.4 Partial Thrust Integrals 


Partial integrals are developed in this subsection, assuming that 
K = Tqq/ 2 (i.e., AK = 0). 

The thrust acceleration vector now assumes the form 
• 

V j = a(t)i^ , where 

Xf = X cos X(t-K) + (y\) sin X(t-K), and 
a(t) = v ex k/( T k “t) or a(t) = a L , 

K = T gq / 2, and 

k = Number of current thrusting phase. 

For constant thrust phases, integration of equation 28 involves 
sine and cosine integrals, which are series expansions, i.e., 
exact closed form integrals do not exist. However, employing 
mean value considerations, thrust integrals are developed that are 
exact for constant acceleration and virtually exact for constant 
thrust. 


Assume that acceleration can be represented as a linear function 
of time, i.e., 

a(t) = A + B (t-K), where A and B are chosen such that 
/ ^ a(t) dt = L and 



t 

f 


o o 


a(s)dsdt = S, where L and S are defined by 


equations 4 and 5. 



1 



Performing the integration in equations 29 and 30 and solving for 
A and B results in 
A = L/Tqq and 

B = 12D/Tg 0 = 3D/2K 3 , where 
D = LK-S. 

• 

Referring to the above equations for Vy, iy, ar.d a(t), the resulting 
thrust acceleration is represented as 

Vj s [A+(B/A) A (t-K)][A. cos A (t-K) + (A/A) sin A f t-k)]* • 
Making a change in the independent variable: 
x = X ( t-K) , 

0 = AK, or A = 0/K, and 

dt = dx/A = (K/0)dx. 

The integration limits are as follows: 

when t=0, x=-0, and 

when t=T G0 , x = 0. 

From the above, it follows that 

dVy = (K/9 )[A+(BK/G )x][_X cos x + (AK/Q) sin x] dx, and 

dRj = (K/0)[/ X dV ]dx. 

-0 

Integration of the above equations results in 

0 0 
V T = / dV, and R T = / dR T , where 

-0 -6 

Vy is velocity change due to thrust and Ry is position 
change due to thrust. 

This results in tv. linear equations of the form: 


*The term in this equation involving the constant, B, is multiplied by 
A and divided by A such that this equation involves terms of the form 
x cos x and x sin x, where x is defined as A ( t - K). 
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and 


R y = SpX + QpX , where 

Lp, Op. Sp, and Q p are thrust integrals developed below. 


( 32 ) 


The following definitions are made: 

A 1 = AK/0 = L/20 
A 2 e A ( K/ 0 ) 2 = LK/29 2 
A 3 e A( K/ 0 ) 3 = LK 2 /20 3 
B 1 = B(K/e) 2 = 3D/2K9 2 
B 2 E B(K/0) 3 = 3D/20 3 
B 3 E B(K/0) 4 = 3DK/20 4 
It follows that 

dVy = C(A ] cos x + B 1 x cos x)X + (A £ sin x + B 2 x sin x)Xjdx, 

X * 

dRy = /•[( A 2 cos s + B 2 s cos s) X + (A 3 sin s + B 3 s sin s) X] dsdx, and 

0 6 
Vy = / dVy, and Ry = / dRy. 


Comparing the above equations with equations 31 and 32, it follows 
that 


0 


0 


L p * A, / cos x dx + B. / x cos x dx 
r -0 '-0 


0 


0 


Op = A ? / sin x dx + B, / x sin x dx 
-0 ^-0 

0 x 0 x 

S p = A o / / cos s dsdx + B 9 / / s cos sds dx 
- 0-0 - 0-0 

Ox 0 X 

Qp = A- / / sin s dsdx + B, / / s sin s dsdx. 
- 0-0 J - 0-0 
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Employing integration by parts, 

/ x cos xdx = cos x + x sin x and 

f x sin xdx = sin x - x cos x. 

The following integrals are readily determined as 

6 6 

/ cos xdx = 2 sin 0, / sin xdx = 0 
-0 -0 

0 0 

/ x cos xdx = 0, / x sin xdx = -2 (0cos0 - sine) 
-0 -0 

0 x 

/ / cos s dsdx = 20 sin 0 
-0 -0 

0 x 

/ / s cos sdsdx = -4 (0 cose - sine) - 26 2 sin 6 
-0 -0 

0 x 

f f sin sdsdx = 2 (6 cos 6 - sine) 

-0 -0 

6 x 

/ / s sin sdsdx = -20 (0 cos 0 - sin 0) 

-0 -e 


From which it follows that 

L p = (L/20) (2 sin9), or 
Lp = L (1/e) sin 0. 

Defining f^ = (1/0) sin 0, then 


The expression for Jp is 

Jp = (3D/20 3 ) [-2 (0 cos 0 - sin 0)] or 

Jp = D |3[(l/0)sin G - cos 0]/O 2 ] } = D[3(f^-cos O)/0 ? ]. 
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Defining fg = 3(f^ - cos 0)/6 , then 
Op = f 2 D. 

3 2 

If it is assumed that sin 6 = 0 - 0 /6 and cos 6 a 1 -0/2, then 

2 

fj = 1 -0 /6 and f ^ = 1 , therefore, it is seen that f^ and f 2 
approach the value of unity as 0 approaches zero. 

The expression for Sp becomes 

Sp = f^LK + D (2f g - 3f ^ ) or 
S p * [LK - D(3-2'f 2 /^*|)3f 1 . and 
Q p = (-LK 2 /3 + DK)f 2 . 


If it is assumed that = 1 in the equation for Sp, then 

Sp = fjS, since D = LK-S. 

This same result is obtained by assuming that a s 2S/Tgg (constant, mean 

value of acceleration for distance) and 

S p = a / G0 / cos X(s-K) dsdt = f,S. 

0 0 

This equation for Sp will be used in the following analysis. 


In summary: 


6 - XT 60 /2 

jp = f 2 D 

f^ = (1/G) sin 0 

S p = f } S 

fo = 3(f^ - cos 0)/O 2 

Qp - (-LT^/12 + DT G0 /2)f 2 , 

Lp - f,L 

since K = Tgg/2. 

0 - LT go /2 - S 
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3.5 Total Thrust Integrals 

Total integrals are developed below, where K = T G q/ 2 + AK(i.e., AK 0). 
The total thrust integrals are now readily obtained, making the 
following definitions: 

6 = XAK, x = X(t - T gq /2) 

cos (x-6) = cos x cos 6 + sin x sin 6 

sin (x-6) - sin x cos 6 - cos x sin 5 

a(x) = A + (BK/0 )x 
0 

Jj = ('70) 2 / a(x) sin (x-6) dx = 0 
-e 

6 x 

Q t = (K/9) 3 f f a ( s) sin (s-6) dsdx 
1 -6 -0 

6 x 

S T = (K/9) 2 f f a(s) cos (s-6) dsdx 
1 -0 -0 

0 

L t = ( K/0 ) / a(x) COS (x-5)dx 
-0 

From inspection of the above equations, it follows that: 

• • 

Oj = Jp cos 6 - (Lp/X) sin 6=0 (since K/0 = 1/X) or 

tan 6 = XJp/Lp = tan (XAK) = X(b/L) (fg/fj ) . 

• • 

If it is assumed that tan (XAK) = XAK* and that ^ 2 ^] ~ ^ ™ the a ^ ove 
equation, then 

AK = D/L = (LT gq /2 - S)/L, or 


*It can be shown that for a single-stage, constant thrust burn, AK 
is closely approximated as 

AK = (T G0 /6)(l-e" V G0 /,V ex)/(l+e‘ V G0^ V ex) , i.e., the upper 
limit for AK is T G q/ 6, and when V GG is small compared to V gx , AKs 0. 

So, 6 = XAK has an upper limit of 0/3. F° r a constant acceleration burn, 
AK=0. For a low thrust OMS burn AK «• 0, and AK converges to the value of 
zero in all cases. So, tan6 = 6 is a valid approximation in all cases. 
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AK = T 60 /2 - S/L, and 
K = T gq /2 ♦ AK = T 60 - S/L, 

which is identical to the current computation of equation 16. 

The expression for becomes 

* 

Qj = Qp cos 6 - (1/x) Sp sin6. 

From the above equation for J-j., it is seen that sin 6 = (XO p /L p ) cos 6, or 
sin 6= (XDf 2 /^ 1 L ) cos 6, from which 

Q t = (-LTqq/12 + DT gq / 2) f 2 cos 6 - (Df^S/f^) cos 6 or 

Q t = (-LTqq/ 12 + DT gq /2 - DS/L) cos 6. (34) 

The above equation can be written as 

Q t « L[-Tg 0 /12 + (AK) 2 ]f 2 cos 6. . 

It follows that 

« • 

Sj = Sp cos 6 + XQp sin 6, where sin 6 = (XJp/L p ) cos 6. 

The second term in the above equation is small compared to the first 
term. If it is assumed that 
Qp ~ SpTeg/6 » snd 

sin 6 5 cos 6, the above equation becomes 

$j = (S p - X6S p T GG /6) cos 6 or 

Sj = S(l-G6/3)fj cos 6 (since XT G q = 2G). (35) 

The same result is obtained by assuming a constant, mean acceleration of 
a(s) = 2 S/T 2 q 

and the above approximation for Q p is not necessary. 
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The integral Ly is simply 
Ly » Lp cos 6 or 
Ly * Lfj cos 6. 

The same equation for Ly is obtained by assuming a constant, 
mean acceleration of a(x) = L/T,,q. 

Equation 34 car >e written as 
Qy c QjFg, where 

Q 1 e ” LT 60 /12 + DT G0 /2 " DS/L * and F 2 = f 2 cos 6 * 

No approximations were made in deriving the above equation other 

than the original assumption of linear acceleration. Similarly, 

if a constant, mean acceleration of 2 S/T^.q is assumed, it can be 

shown that Qy = Qq F 2’ where Q() = ”$ T G(/ 6 * DS/L. The significance 

of the above equations for Qy is that Qq and Q-j are independent 

of the value of X, and the function F 2 factors out of Qq and Q-j. 

Simulations have shown that, in general, Q q is a fairly crude 

approximation for the first order integral, Q-SK. Qj is a sig- 

nifant improvement over Qq, however, an error of approximately 5% 

can exist, e.g., during the first PEG phase of Baseline Reference 

Mission 3A.* This reflects an error of approximately 5# in X_ 

and 6, since, from equation 27, it is seen that 

* • 

A. = (£gn " S-y?L)/Qj * i-e., X. is inversely proportional to Q^ . 

*In these simulations, Qy was formed as a multi-stage integral, 
analogous to the multi-itage value of Q-SK. 
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Similar to the above considerations, a second order acceleration 
profile can be assumed, i.e., 

a(t) = A + (B/X)X(t-K) + (C/X 2 )X 2 (t-K) 2 
and the expression for Qy can be expressed as Qy s Q 2 F 2 , where 
Qg is independent of the value of X, and is a very good approximation 
for Q-SK. Closed form integrals exist for any order of a ( t) , 
however, the expression Q n becomes complicated and the integral 
Q-SK is simpler to implement. By the process of induction, the 
integral Qy can be expressed as 

Qy = (Q-SK) F 2 > where F 2 * f 2 cos 6. 

The above equations are summarized as: 

9 = i T gq /2 
f-j = (1/0) sin 0 
f 2 = 3(f 1 - cos 0)/6 2 
K = T G0 - S/L 
S * X(K-T go /2) 

» f^ cos 6 
F 2 = f 2 cos 6 

Fj = F^ (1 -06/3) (F^om equation 35) 



q t = f 2 (q-sk) 

where L, S, and Q are the first order, multi-stage integrals 
defined in Section 3.1. 
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The results of the above analysis can be summarized as follows: 

a) Defining AK « K-Tqq/2, where K is the ci „».unt, reference 
time aoout which the solution Is expanded, It was shown 
that the resulting value of AK Is the same whether a mean, 
linear acceleration profile is assumed or the actual 
acceleration profile Is assumed. It can also be shown that 
(assuming linear acceleration or actual acceleration) the 
multi-stage expression for AK Is the same as the single- 
stage expression. It can be concluded that the value of AK 

is directly proportional to the mean slope of the acceleration 
profile (multi-stage or single-stage), since AK = D/L and 
the expression for the mean slope is B - 12 D/Tqq. 

b) The total integrals Ly and Sy are obtained by assuming 
a constant, mean acceleration (a^L/T^g for deriving Ly, 
and a = 2 S/T£q for deriving Sy). 

c) The integral Qy is a higher order i teg^al and approaches 

the value of (Q-SK^ as the order of the assumed acceleration 
profile increases. 

d) Since a constant acceleration can be assumed in deriving 
the integrals Ly and Sy, it can be concluded that the 
factors F.| and F^, as defined above, apply to an entire 
maneuver, i.e., do not have to be derived in a multi-stage 
manner. Since F^ * 1, it will be assumed that the above 
conclusion also applies to this factor. 
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e) It was shown in Section 3.1 that the expression for the 
integral assuming constant acceleration, is the same 
as that when assuming constant thrust. This simplifies 
Implementation of the thrust integrals, since two different 
expressions are unnecesssary. 

f) The higher order thrust integrals in the current PEG 
program can be eliminated by introduction of the factors 

, and F^ . 
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3.6 Implementation of Thrust Integrals 

The proposed thrust integrals are presented in the flow chart of 
Figure 3.6-1. For purposes of comparison, the current thrust 
integrals are presented in Figure 3.6-2. The proposed integrals 
require about one half the code required by the current integrals 
and approximately 30 data words (assuming double precision) are saved by 
eliminating the variables S^, J., Q^, P^» and H.. (i = 1, 3). 

Rationale for implementing the integrals as in Figure 3.6-1 was 
discussed in Section 3.1 and additional discussion is included 
in this section. The overall method of implementing the thrust 
integrals is presented in this section. 


Currently the desired unit thrust vector is computed in the 
Guidance and Control (G&C) steering interface routine, not in the 
PEG routine. The parameters and T^ = Tq + K are sent to the G&C 

steering interface from PEG (where Tq is the guidance time when K was 
computed). The unit thrust vector is then computed as 


Vp = unit [X + (t-T^)Xj , 
which is identical to 


(36) 


Vp = unit [X - KX + (t-Tg)Xj, where the (37) 

value of t-Tg ranges from near zero to the value of the guidance 

cycle time (e.g., 2 seconds). Equation 21 of Section 3.2 can be 

written in an identical form to equation 37 by redefining the parameter 

T^, which is currently the last computation in the thrust integrals. 

Equation 21 can be used to define ip as ip = unit [X cos uz + (yw) sin wz]? 

where z = t-Tg-K. The above equation can be written as 
* At this point, the parameter X is replaced by w, since in general 
(a X. 
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Xf * unit [X. + X (jjj tan wz)] SIGN (cos wz). 

If the angle wz is limited to be less than it/2, SIGN (cos wz) * +1, 
and since w(t-Tg) is a small angle, the above equation can be closely 
approximated as 

i.f = unit [X - K(1/wk)(tan wk)X + (t-Tg)X]. 

The above equation can be written as 

i^ = unit [X. - KpX + (t-Tg)Xj, where 
K p = F 4 K, and F 4 2 (1/wK) tan wK). 

Defining as + K p , the above equation becomes 

= unit [X + (t-T x )X], 

which is identical to equation 36. Implementation of the thrust 
integrals in the above manner has no impact on the current G&C 
steering interface equations, and only one scalar equation is 
modified in the PEG thrust integrals. 

As discussed in Section 3.2, using linear angle steering (i.e., 

» 

to = X) produces a MECO weight loss for Baseline Reference Mission 
3A. Experience has shown that linear tangent steering is more 
optimum for Mission 3A. The steering is linear tangent if 

w ~ 0, e.g., w = .00001. This capability is implemented by initial - 

izing w = .00001 (in the PEG initialization block) and making the 
following test in the PEG turning rate block. 

If w < .00001 or n=3, w = .00001, where n is the number of phases, and 

Mission 3A is the only Shuttle maneuver with as many as 3 phases. The 

above test without "or n * 3" is equivalent to: If w = 0, w > 0. 
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This test is recommended because the parameter fj is defined as 
fj= |sine, where 0 = wTqq/ 2, i.e., division by zero is avoided. 

This results in linear tangent steering and linear tangent prediction. 
Simulations have shown that this is more optimum for Mission 3A 
(+70 pounds MECO weight) than the current linear tangent steering 
and linear sine prediction. 


Experience has shown that when the thrust integrals are computed in 
single precision, they become unstable during the last 40 seconds of 
an OMS burn. This is probably due to instability in the 0 and Q 
integrals. These integrals can be written as 

■> ■ ■'(‘--“oW and 

Q = x(S-i a T^q) , where a is the value of current 
acceleration. The coefficients of x, in the above equations, approach 
the value of zero as Tg Q approaches zero. Single precision 
probably produces randomness in these coefficients. If acceleration 
is expanded as a quadratic, the above equations are approximated as 

j *K t go< ,+ ^¥> 

O'IVgo * 1 + k;¥»- 

The value of T^q/x approaches zero as T^q approaches zero, and the 
expressions for J and Q approach J = ^ a Q TgQ and Q = ^ a o^G0* 

If an average acceleration of V^q/T^q is assumed the values of J and 
Q can be approximated in a virtually exact manner during the last 
40 seconds of an OMS burn, i.e.. 
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0 = 2 V GO T GO and 

Q = ^ VqqTqq (at this time AK » 0). 

Single precision is adequate if the equations for J and Q are expressed 
as above, and this result is accomplished if the constant acceleration 
thrust integrals are used when T^q is less than 40 seconds (for an 
OMS burn). This capability is implemented in the thrust integrals 
of Figure 3.6-1. 

The equations for the total thrust integrals, as summarized in 
Section 3.5, produce more accuracy than is necessary in practice. 

For example, assumptions used for Saturn V/Apollo guidance were 
F^ = Fg = F 3 = 1 . However, in cases where accurate prediction is 
desirable (e.g., throttling is involved or deorbit targeting) the 
F^ term is important. The angle 6 is at most about 4 degrees 
during shuttle second stage burn, therefore it can be assumed 
that cos 6=1. The quantity 66/3 is at most about .01, and can 
be neglected*. As observed in Section 3.4, = 1 for third order 

assumptions on sine and cosine expansions. Also, this term has 
no direct effect on velocity prediction, therefore, it can be 
assumed that f ^ = 1 . 

Experience has shown that the most stable manner of implementing the 
f-j term is to feed back information from the previous guidance 
cycle. Currently the quantity "* s f ec * back. The ^ erm 

can be eliminated by initializing f-j = 1 and defining S * f^S. 

♦However, if throttle control is used, it may be desirable not to 
neglect this term. 
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Rqq is formed without the Rg j AS term as R^q = R^, where R^ is 
defined as in Section 3.3. The above form is somewhat more 
accurate than current implementation of the Rgj^<. term. S is 
computed in the range-to-go block and f^, in the predictor block. 


The V. 7 HRUST com P utatlon and variable name is eliminated in the predictor 
block by replacing this quantity with f Ao’ and JLjhrusT re P^ aced 
with Rqq, where f-j has the value of the current guidance pass. 


Currently, the thrust integral Q-SK is used at six places in the 
PEG equations, and recomputed each time. This is avoided by 
replacing Q-SK with Q = Q-SK. 


The test, if Q-SK f 0, can be eliminated in the turning rate block 
since Q-SK can not be zero unless Tqq = V G q = 0, and if this is the 
case the guidance routine has already been aborted. 


The <j>|^ equations are simplified by elimination of the equations for 
recomputing \ and ^GO’ since the main purpose for this is for the 
computation of and this computation is eliminated. However, 

it is desirable to limit the value of f-j, and this is accomplished 
by maintaining the test. 
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The above discussion is summarized below and detailed implementation 
is discussed. Modifications required in each PEG block are presented. 
An estimation of code saved in each block is included. 

Initialization Block 
Eliminate: = 0 

Add: fj = 1 and to = .00001 

(code saved: =0 words) 

Integrals Block 

Replace the equations of Figure 3.6-2 with the equations 
of Figure 3.6-1 . 

(code saved: 100 words) 

Turning Rate and Range-To-Go-Blocks 
Replace Q-SK with Q. 

(code saved: 30 words) 

Turning Rate Block 
Eliminate the following equations: 

If Q-SK t 0, Then X = (R^ - SX)/(Q-SK) 

, Else X_ = 0 

X = lii 

If K ^mr * = 

• • • 

, X. = X unit (X) 

, ^ = SX + (Q-SK) X 
Replace the above equations with: 
i « (Bgo - SX)/Q 

w = lii 
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If u> * <t> MAX / K 

If a) < .00001 or n s 3, w * .00001 

(code saved: 50 words) 

Range-To-Go Block 

Define: ^ = Rq - (R + V T G0 + Rg rav ). 

Replace R^q = Rq N + Rg IAS with 

^G0 = ^GN and S " f l S * 

(code saved: -10 words) 

Predictor Block 

Eliminate: 

-THRUST = f- L " 2 

—THRUST = - \ ^ 2 (P- 2 Q K + SK 2 )]X + (Q-SK)i 

^BIAS = ^GO ‘ -THRUST 
Replace with: 

f l = (2/wTgQ) sin (coT gq /2) 

Replace Rj HRUS T with ^GO and -THRUST with f l ^GO* 

(code saved: 50 words) 

The total estimated code saved from the above PEG modifications 
is 220 words. This amounts to 10 percent of the current PEG 
program. 
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For OMS maneuvers, involving no position constraints, the thrust 
integrals could be eliminated completely by sending the total 
acceleration (thrust and gravity) to the average G routine, returning 
with predicted velocity and position. An OMS (one stage) version of 
PEG would be extremely simple implemented in this manner. In 
general, analytical thrust integrals could be eliminated by use of 
numerical integration, however, implementation is more difficult if 
position constraints are involved. 

The most accurate and optimum and possibly the simplest manner to 
implement PEG is to update Rg 0 in the PEG update block as 
^60 = ^GO A V T G0 - Atg/2] and recompute the parameter in the PEG 
correction block as Rg 0 = + Bq “ Bp* Then, the thrust acceleration 

vector, (F/m)i^, where i^ is a unit vector, is sent to the average G 
routine? along with the gravity vector, and values of predicted 
velocity and position are obtained. This eliminates the requirement 
for gravity integrals and higher order thrust integrals, since 
these quantities are implicit in the resulting values of ho and 
^GO* The entire PEG predictor block would consist of defining a 
value of the thrust acceleration vector and the integration time step 
and calling the average G routine? This would entail multi-stage ' 
logic and logic for defining the integration time step. If PEG 
were implemented in this manner, the term involving sin 0 defined 
in Subsection 3.2, could be incorporated into the thrust acceleration 
vector, and the PEG program would be equivalent to an accurate 
onboard calculus of variations trajectory optimization program. 

*or some other numerical integration routine. 
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An alternate consideration is to compute R^ 0 as defined above and 

maintain the current PEG predictor method. However, the variable 

names and equations for V and R can be eliminated and the 

-Kjrav -=qrav 

equations for predicted velocity (V^) and position (R^) can be 
written as 


V Sci ♦ f l *» ♦ * - *el and 


= ^c2 + ^cl " - + f l 4oW 6, 


In addition, the variable names hz and ^c2 can be eliminated by 
replacement with Vp and R^. This method of implementation would 
result in a savings of approximately 50 words, in addition to the 
220 words discussed earlier in this section. 
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4.0 CONCLUSIONS 

The new PEG thrust integral formulation and method of implementation 

derived and presented in this note have the following advantages: 

a) Economy - Preliminary estimates' indicate a computer storage 
savings of 220 words, approximately 10 percent of the current 
PEG program. 

b) Ease of implementation - The implementation method proposed in 
this note does not represent a departure from the current PEG 
implementation approach, rather a simplification of current 
equations. The new equations are implemented with very minimum 
effort. 

c) Ease of verification - The implementation approach (logic, 
sequence of computations, etc.) requires no verification, since it 
is analogous to current implementation, i.e., stability and con- 
vergence are not affected. The first order integrals of the 
current PEG equations are maintained, therefore, require no verifi- 
cation. The simplified higher order integrals are verified 
analytically. The primary motive for verification is to demonstrate 
the additional flexibility. 

d) Flexibility - The proposed equations represent a generalized 
guidance law, i.e., different guidance laws are applied by 
changing the value of one scalar parameter. This offers opti- 
mality and accuracy of prediction for various types of maneuvers. 

e) Accuracy - The new equations produce more accuracy than current 
equations if large central turning angles are Involved (e.g., 
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greater than 90 degrees), since sine and cosine approximations 
are eliminated. 

f) There are no known disadvantages to the proposed equations. 

Altetnate implementation methods, which could result in even more 
computer storage savings, are noted for possible consideration. 
However, the alternate methods represent a larger deviation from the 
current PEG equations than the primary method presented. 

The proposed equations were verified in simulations of Baseline 
Reference Missions 1 and 3A. Slightly improved performance was 
demonstrated for high performance missions such as Mission 3A 
(i.e., +70 pounds weight at MECO). 

It is recommended that these equations be incorporated into the 
Shuttle powered explicit guidance software. 
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